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Annihilation of Dark Matter usually produces together with gamma rays comparable amounts of 
electrons and positrons. The e^e~ gyrating in the galactic magnetic field then produce secondary 
synchrotron radiation which thus provides an indirect mean to constrain the DM signal itself. To 
this purpose, we calculate the radio emission from the galactic halo as well as from its expected 
substructures and we then compare it with the measured diffuse radio background. We employ a 
multi-frequency approach using data in the relevant frequency range 100 MHz-100 GHz, as well as 
the WMAP Haze data at 23 GHz. The derived constraints are of the order {aAv) = 10~'^*cm^s~^ 
for a DM mass ~ 100 GeV sensibly depending however on the astrophysical uncertainties, in 
particular on the assumption on the galactic magnetic field model. The signal from single bright 
clumps is instead largely attenuated by diffusion effects and offers only poor detection perspectives. 

PACS numbers: 95.35.+d, 95.85.Bh, 98.70.Vc 



I. INTRODUCTION 

Cosmology and Astrophysics provide nowadays a com- 

Silhng evidence of the existence of Dark Matter (DM) 
,0]. Nevertheless, its nature still remains elusive, and 
Dark Matter constituents have escaped a direct detec- 
tion in laboratory so far. Promising candidates are DM 
particles produced in thermal equilibrium in the early 
universe, the so-called Weakly Interacting Massive Par- 
ticles (WIMPs). Theoretically, models of WIMPs natu- 
rally arise, for example, in SUSY as the Lightest Super- 
symmetric Particle or as the Lightest Kaluza-Klein Par- 
ticle in the framework of extra-dimensions. These candi- 
dates are self-conjugate and can thus annihilate in cou- 
ples to produce as final states: neutrinos, photons, elec- 
trons, light nuclei (as wells as their antiparticles), etc., 
which can in principle be detected. 

Among the indirect DM detection channels, gamma- 
ray emission represents one of the most promising oppor- 
tunity due to the very low attenuation in the interstellar 
medium, and to its high detection efficiency. See for ex- 
ample Ref.s 0,1,13 for a review of this extensively stud- 
ied issue. The expected neutrino detection rates are gen- 
erally low although forthcoming km'^ detectors offer some 
promising prospect 0] . Finally, positrons and protons 
strongly interact with gas, radiation and magnetic field 
in the galaxy and thus the expected signal sensibly de- 
pends on the assumed propagation model 0, [8|, [Sj] • How- 
ever, during the process of thermalization in the galactic 
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medium the high energy e+ and e~ release secondary 
low energy radiation, in particular in the radio and X- 
ray band, that, hence, can represent a chance to look for 
DM annihilation. Furthermore, while the astrophysical 
uncertainties affecting this signal are similar to the case 
of direct e"*", e~ detection, the sensitivities are quite dif- 
ferent, and, in particular, the radio band allows for the 
discrimination of tiny signals even in a background many 
order of magnitudes more intense. 

Indirect detection of DM annihilation through sec- 
ondary photons has received recently an increasing at- 
tention, exploring the expected signature in X-rays [Tol . 
[Ill, [l^ ■ at radio wavelengths jl3, 14, J5|, , or both 
[ITI. llSl. [Toj . In the following we will focus our anal- 
ysis on the radio signal expected from the Milky Way 
(MW) halo and its substructures. It is worth noticing 
that the halo signal has been recently discussed in Ref.s 
[13, mi. Hi in connection to the WMAP Haze, which has 
been interpreted as a signal from DM annihilation. In 
this concern we will take in the following a more con- 
servative approach, by assuming that the current radio 
observations are entirely astrophysical in origin, and thus 
deriving constraints on the possible DM signal. The main 
point will be the use of further radio observations besides 
the WMAP ones, in the wide frequency range 100 MHz- 
100 GHz, and a comparison of the achievable bounds. 
Furthermore, the model dependence of these constraints 
on the assumed astrophysical inputs will be analyzed. 
We will also discuss the detection perspectives of the sig- 
nal coming from the brightest DM substructures in the 
forthcoming radio surveys. 

The paper is organized as follows: in section |TT] we will 
discuss the astrophysical inputs required to derive the 
DM signal such as the structure of the magnetic field, 
the DM spatial distribution and the radio data employed 
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to derive the constraints. In section IIIII we describe in 
detail the processes producing the DM radio signal ei- 
ther when it is originated from the halo or from the sub- 
structures. In section IV we present and discuss our con- 
straints, while in section V we analyze the detection sen- 
sitivity to the signal coming from the single DM clump. 
In section VI we give our conclusions and remarks. 



II. ASTROPHYSICAL INPUTS 

A. Dark matter distribution 

Our knowledge of the DM spatial distribution on galac- 
tic and subgalactic scales has greatly improved thanks 
to recent high resolution zoomed N-body simulations 
[H, H^j H^, l26l | . These simulations indicate that for 
the radial profile of the galactic halo the usual Navarro- 
Frank- White (NFW) distribution [13 



pir) = 
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still works as a good approximation over all the resolved 
scales. We will thus use this profile in the following. 
Note, anyway, that this choice is quite conservative with 
resp ect to other proposed profiles like the Moore profile 
[2^ . which exhibits an internal cusp oc r"^-^ that would 
give in principle a divergent DM annihilation signal from 
the center of the halo. Observationally the situation is 
more uncertain: Baryons generally dominate the gravi- 
tational potential in the inner kpc's and fitting the data 
thus requires to model both the baryon and DM com- 
ponent at the same time. The NFW profile is in fair 
agreement with the observed Milky Way rotation curve 
[29( 1 ■ although, depending on the employed model, it is 
possible to find an agreement for many different DM pro- 
files (see also Q and references therein). We emphasize, 
however, that the various profiles differ mainly in the 
halo center (for r < 1 kpc) where the uncertainties, both 
in numerical simulations and from astrophysical observa- 
tions are maximal. Thus, our analysis which explicitly 
excludes the galactic center, does not crucially depend 
on the choice of the profile. 

A problem related to the profile of Eq[T]is that the mass 
enclosed within the radius r is logarithmically divergent. 
A regularization procedure is thus required to define the 
halo mass. Following the usual conventions we define 
the mass of the halo as the mass contained within the 
virial radius rvir, defined as the radius within which the 
mean density of the halo is S^ii = 200 times the mean 
critical cosmological density pcr which, for a standard 
cosmological model (fim — 0.3,17a — 0.7 [ij) is equal to 
Per — 5 X 10~^GeVc~^ cm~'^. The parameters describing 
the halo arc then determined imposing the DM density 
to be equal to ps = 0.365 GeV cm~^ near the Solar 
System, at a galactocentric distance of i?s = 8.5 kpc. 



Simulations, however, predict a DM distribution sum 
of a smooth halo component, and of an additional clumpy 
one with total masses roughly of the same order of magni- 
tude. Hereafter we will assume for the mass of the Milky 
Way Mmw = Mi + A^ci = 2 x IO^^A/q, where A'h and 
Aid denote the total mass contained in the host galactic 
halo and in the substructures (subhalocs) distribution, 
respectively. The relative normalization is fixed by im- 
posing that subhaloes in the range lO^M© , lO^^M© have 
a total mass amounting to 10% of Mmw [111- Current 
numerical simulations can resolve clumps with a mini- 
mum mass scale of ~ 10^ Af©. However, for WIMP par- 
ticles, clumps down to a mass of 10~^AfQ are expected 
[sol , [sij . We will thus consider a clump mass range be- 
tween 10-^ Mq and 10^^ Mq. 

Finally, to fully characterize the subhalo population 
we will assume a mass distribution oc m~^^ and that they 
are spatially distributed following the NFW profile of the 
main halo. The mass spectrum number density of sub- 
haloes, in galactocentric coordinates r, is thus given by 



(mci,r) = A — 

drrici \Mc\J V^h 
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where A is a dimensional normalization constant. The 
above expression assumes some approximations: for ex- 
ample, a more realistic clump distribution should take 
into account tidal disruption of clumps near the galac- 
tic center. Numerical simulations suggest also that the 
radial distribution could be somewhat anti-biased with 
respect to the host halo profile. However, with our con- 
servative assumptions the DM annihilation signal is dom- 
inated by the host halo emission within up to 20° — 30° 
from the galactic center so that the details of the sub- 
dominant signal from the clumps have just a slight influ- 
ence on the final results. Recent results also show that 
mass distribution seems to converge to 'n\^''^ rather than 
TT^cX^ [25j . This would also produce only a minor change 
in the following results. 

Following the previous assumptions the total mass in 
DM clumps of mass between mi and m2 results to be 
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where Ch = 7'vir/ni denotes the host halo concentration; 
while their number is 
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Imposing the normalization condition 

M{10HIq,10^°Mq) = 10%Mmw, we finally get 
for the mass due to the entire clumps distribution: 



Radio Pkv at 23 GHz 



Mcl = M(1O-^M0, lO^^Mo) - 53.3%AfMW 



while for the number of these clumps we obtain 



NilQ-^MQ, 10^°Mq) ~ 2.90 x 10 
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(5) 



(6) 



Finally by using the previous constraints one can fix the 
values of free parameters rh, Ph and A, hence obtaining 
Th = 14.0 kpc, that corresponds to a halo concentration 
of Ch = 14.4, Ph = 0.572 GeV c'^ cir-^ and A = 1.16 x 
10-19 kpc-3 M-\ 

A further piece of information is required to derive 
the annihilation signal from the clumps, namely how the 
DM is distributed inside the clumps themselves. We 
will assume that each clump follows a NFW profile as 
the main halo with and pci replacing the correspond- 
ing quantities of Eq[TJ However, for a full character- 
ization of a clump, further information on its concen- 
tration Cci is required. Unluckily, numerical simulations 
are not completely helpful in this case, since we require 
information about the structure of clumps with masses 
down to 1O-^M0, far below the current numerical res- 
olution. Analytical models are thus required. In the 
current cosmological scenario [l| structures formed hier- 
archically, via gravitational collapse, with smaller ones 
forming first. Thus, naively, since the smallest clumps 
formed when the universe was denser, a reasonable ex- 
pectation is Cci (X {1 + Zf), where z/ is the clump for- 
mation redshift. Following the model of ref. |32| we 

will thus assume Cd = ci (^jf^^ with ci ~ 38.2 and 

a = 0.0392. With this concentration the integrated DM 
annihilation signal from all the substructures dominates 
over the smooth halo component only at about 30° from 
the galactic center (see section IIIip . so that the con- 
straints on the DM signal do not crucially depend on the 
unresolved clumps signal, coming basically only from the 
smooth halo component. However, given the large un- 
certainties in the models, larger contributions from the 
unresolved population of clumps are in principle possible 
considering a different parametrization of the concentra- 
tion (see for example the various models considered in 
0, IS]). We will not investigate further this possibility 
here. An enhancement of the clumps signal is also possi- 
ble considering different choices of the clump profile other 
than the NFW: Differently from the case of the halo, in 
fact, the clump signal depends sensibly from the chosen 
profile and a Moore profile or an Isothermal profile can 
in principle enhance the signal of several orders of mag- 
nitude. Also in this case we choose to quote conservative 
constraints and we will not consider these possibilities 
further. 
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FIG. 1: Sky map of galactic foregrounds at the frequency of 
23 GHz (top), and of the residual map showing the WMAP 
Haze (bottom). 



B. Galactic Magnetic Field 

The MW magnetic field is still quite uncertain espe- 
cially near the galactic center. The overall structure 
is generally believed to follow the spiral pattern of the 
galaxy itself with a normalization of about ~ IpG near 
the solar system. Eventually, a toroidal or a dipole com- 
ponent is considered in some model. 

We will consider in the following the Tinyakov and 
Tkachev model (TT) [sj] which is a fair representative of 
the available descriptions of MW magnetic field. Within 
this model , the field shows the typical spiral pattern, an 
exponential decrease along the z axis and a 1/ R behavior 
in the galactic plane. The field intensity in the inner kilo- 
parsecs is constant at about 7fj,G. We will use the slightly 
modified parametrization of this model as described in 
(3^ . Higher normalizations are in principle possible con- 
sidering more complex structures as for example a dipole 
or a toroidal component [3^ . Indeed some recent anal- 
yses [13, nil including new available data seems to favor 
the presence of these further structures. We will thus 
consider as possible also an "high normalization model" 
that we simply parameterize as a constant 10 fiG field. 
This choice is also motivated by a comparison with the 
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results of [20|, |2l| where the same magnetic field is used. 

Further, beside the regular component, the galactic 
magnetic field presents a turbulent random component. 
The r.m.s. intensity of this component is generally ex- 
pected of the same order of magnitude of the regular one, 
but both its spatial distribution and its spectrum are 
poorly known, thus here we neglect its effects. Naively, 
this random component is expected to affect the syn- 
chrotron maps that we will show in the following pro- 
ducing a blurring of the otherwise regular pattern. Also, 
the random component contributes to increase the overall 
normalization of the field. Thus without this component 
the synchrotron signal is slightly underestimated so that 
we can regard this choice as conservative. 



C. Radio Data 

In the following we will derive constraints on the DM 
emission comparing the expected diffuse emission from 
the smooth halo and the unresolved population of clumps 
with all sky observation in the radio band. In the fre- 
quency range 100 MHz-100 GHz where the DM syn- 
chrotron signal is expected, various astrophysical pro- 
cesses contribute to the observed diffuse emission. Com- 
peting synchrotron emission is given by Cosmic Ray elec- 
trons accelerated in supcrnovae shocks dominating the 
radio sky up to 10 GHz. At higher frequencies the Cos- 
mic Microwave Background (CMB) and its anisotropics 
represent the main signal. However, thanks to the very 
sensitive multi-frequency survey by the WMAP satel- 
lite, this signal (which represents thus a background for 
DM searches) can be modeled in a detailed way and can 
thus be removed from the observed radio galactic emis- 
sion [3^. Other processes contributing in the 10-100 
GHz range are given by thermal bremsstrahlung (free- 
free emission) of electrons on the galactic ionized gas, and 
emission by small grains of vibrating or spinning dust. 

In the following our approach will be to compare the 
DM signal with the observed radio emission where only 
the CMB is modeled and removed. For this purpose wc 
use the code described in [i^ where most of the radio sur- 
vey observations in the range 10 MHz-100 GHz are col- 
lected and a scheme to derive interpolated, CMB cleaned 
sky maps at any frequency in this range is described. 

A more aggressive approach would be of course to try 
to model and subtract also the remaining emissions (syn- 
chrotron, free-free, dust) in order to compare the ex- 
pected DM signal with the residual radio map. This 
is indeed the approach followed in [U, [4^ where resid- 
ual maps at the 5 WMAP frequencies are derived using 
spatial templates for the various expected astrophysical 
components. The residual maps then exhibit the fea- 
ture called the WMAP Haze, which has been indeed in- 
terpreted as radio emission related to DM annihilation 
[UEl]. However, the fit procedure used for the Haze ex- 
traction is crucial, and using more degrees of freedom to 
model the foregrounds as performed by the WMAP team 



[44| fails in finding the feature. We will anyway show in 
the following for comparison the constraints derived us- 
ing the Haze residual map at the WMAP frequency of 23 
GHz [4^ . A map of the Haze and of the 1 GHz emission 
is shown in FiglT] We will sec however that within our 
conservative approach comparable or better constraints 
can be obtained thanks to the use of multi-frequency in- 
formation. For a given DM mass, in fact, 23 GHz is 
generally not the best frequency to use and better con- 
straints are instead obtained using observations at lower 
frequencies even without further foreground modeling. 

Definitely, a detailed foreground modeling at all radio 
frequencies would clearly give much stronger constraints 
on the DM signal and/or eventually confirm the DM na- 
ture of the WMAP haze. To this purpose consistent 
progress will be achieved in the next years with the new 
high quality data coming from the PLANCK mission and 
from low frequency arrays like LOFAR and SKA. 



III. DM SYNCHROTRON SIGNAL 
A. Particle Physics 

In a standard scenario where WIMPs experience a non 
exotic thermal history, a typical mass range for these 
particles is 50 GeV < < 1 TeV, while a simple esti- 
mate for their (thermally averaged) annihilation cross 
section yields (ctav) = 3 x 10~^'^cm^s~^/ ilcdmh'^ 01, giv- 
ing (aAv) « 3 X 10~^^cm^s~^ for ricdmh"^ ~ 0.1 as result- 
ing from the latest WMAP measurements [l|. However, 
this naive relation can fail badly if, for example, coanni- 
hilations play a role in the WIMP thermalization process 
[4^, and a much wider range of cross sections should be 
considered viable. In this work we consider values of 
from about 10 GeV to about 1 TeV, and (iJAf) in the 
range (10-26-10-2i)cm3s-i 

The e~^e~ annihilation spectrum for a given super- 
symmetric WIMP candidate can be calculated for exam- 
ple with the DarkSUSY package [47[. However, the final 
spectrum has only a weak dependence on the exact anni- 
hilation process with the channels XX ~^ ZZ, W~^W~ , qq 
giving basically degenerate spectra. For leptonic chan- 
nels like the t~^t~ decaying mode the spectrum differs 
significantly, although this channel has generally a quite 
low branching ratio. For simplicity we will assume here- 
after full decay into qq channel, hence e~ (e"*") will be 
emitted by decaying muons (anti-muons) produced in pi- 



ous decays. In this framework, the resulting e"*", e spec- 
trum can be written as a convolution, namely 



dE^W^{E^)^^^^{E^) (7) 
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with ^ = {rn^/rriTr) , where 



diV, 



(m) 



hat''"' 1 



5 i f E,V 2 f E,'^^ 



6 2 V£^- 



W^{E^) 



E^ ml - ml 



1 15 / TO^C 



m^c^ 16 



3 V^„ 



,(8) 
(9) 
(.10) 



In particular, Eq. ([8]) is the electron (positron) spectrum 
produced in the muon (anti-muon) decay /i~ — > e~v^j,De 
[ji^ — > e'*'f'^;^e)- Eq- ® stands for the (/^^) spec- 
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decay process, and, 



trum from tt" 

finally, Eq. (fTO)) provides a reasonable analytical approx- 
imation of the spectrum of pions from qq hadronization 
[4^. It is worth noticing that to be more accurate Eq. 
([T0| should be substituted by a numerical calculation, 
which however results not necessary for the aim of the 
present paper as discussed in the following. 

In this approximation the final electron (positron) 
spectrum can be cast in a simple polynomial form of the 
ratio i^e/m^c^: 



dE, 



„2 Z^'^J 



E,, 



(11) 



where J = { — I, —5, 0, ^,2, 3} and the coefficients aj are 
listed in Tablell 

The main advantage of using the above analytical ap- 
proximation instead of a more accurate numerical input 
is that, as will be clear in the next section, most of the 
observables for the radio emission will be expressed in an 
analytical form as well. This, in turn, is of help for a bet- 
ter understanding of the physical results. Nevertheless, 
the difference with the complete numerical calculation 
turns to be small, arising only for quite low electron en- 
ergies, and thus for very low radio frequencies. At low 
energies, in fact, the analytical form has an asymptotic 
behavior oc E~^'^ while the numerical spectrum has a 
turn down. From a comparison with the numerical out- 
put from DarkSUSY for a 100 GeV WIMP with 100% 
branching ratio into bb the analytical form is a fair ap- 
proximation until Ee ~ 1 GeV, which for a magnetic field 
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TABLE I: ai values 
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0.456 
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-6.77 
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-0.185 



_B ~ 10 /iG translates into a minimum valid frequency 
V — 10-100 MHz, thus below the frequency window we 
are going to explore (see Eqi20l below) . 



B. Electrons equilibrium distribution 

Dark matter annihilation injects electrons in the 
galaxy at the constant rate 



1 (p{r) 
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On the other hand, the injected electrons loose energy in- 
teracting with the interstellar medium and diffuse away 
from the production site. In the limit in which convec- 
tion and reacceleration phenomena can be neglected, the 
evolution of the e+e" fluid is described by the following 
diffusion-loss equation [il, HO, [5l| 
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where due/dE^ stands for the number density of e+, 
per unit energy, K{Ee,r) is the diffusion constant, and 
b{Ee,r) represents the energy loss rate. The diffusion 
length of electrons is generally of the order of a kpc (see 
section |V| thus for the diffuse signal generated all over 
the galaxy, and thus over many kpc's, spatial diffusion 
can be neglected. This is not the case for the signal 
coming from a single clump for which the emitting region 
is much smaller than a kpc. We will further analyze this 
point in section |Vl By neglecting diffusion, the steady 
state solution can be expressed as 



dE', Q{E',,r), 



(14) 



where r = Ee/b{Ee, r) is the cooling time, resulting from 
the sum of several energy loss processes that affect elec- 
trons. In the following we will consider synchrotron emis- 
sion and Inverse Compton Scattering (ICS) off the back- 
ground photons (CMB and starlight) only, which are the 
faster processes and thus the ones really driving the elec- 
trons equilibrium. Other processes, like synchrotron self 
absorption, ICS off the synchrotron photons, e~^e^ an- 
nihilation, Coulomb scattering over the galactic gas and 
bremsstrahlung are generally slower. They can become 
relevant for extremely intense magnetic field, possibly 
present in the inner parsecs of the galaxy p^ . and thus 
will be neglected in this analysis. 

For synchrotron emission the energy loss is given by 
(for ex. see [13) b^yniEe) = 4/3 caTl^ P^Ub with Ub = 
B^/2fiQ the magnetic energy density so that the time 
scale of the energy loss is: 
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X (kpc) 3x10 cm"^ s ^ , while the relevant clumps parameters are 

BBBBBB shown in table In] 



FIG. 2: Projections of the galaxy in the xy and xz planes 
showing the fractional synchrotron contribution to the e^e~ 
total energy losses for TT model [11] of GMF and Galprop 
model [H^l of ISRF. The synchrotron losses contribute up most 
to 20 % reaching its maximum at the center of the magnetic 
arms. In the remaining regions, included the galactic center, 
ICS is dominating. 



syn 

Similarly, for Inverse Compton emission the energy 
loss is given by 6ics(-E'e) = 4/3 c(TT7^/?^f^rad- The rel- 
evant radiation background for ICS is given by an extra- 
galactic uniform contribution consisting of the CMB with 
UcwB = %TT^{kTf/lb{hcf K. 0.26 cV/cm^, the opti- 
cal/infrared extragalactic background and the analogous 
spatially varying galactic contribution, the Interstellar 
Radiation Field (ISRF). For the latter we use as tem- 
plate the Galprop distribution model [11] which reduces 
to the extragalactic one at high galactocentric distances. 
In this model, the ISRF intensity near the solar posi- 
tion is about 5 eV/cm^, and reaches values as large as 
50 eV/cm^ in the inner kpc's. With this model the ICS 
is always the the dominant energy loss process, also near 
the galactic center (see Fig (2). We thus have 
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with Tj^cg = 9.82 X IQi^s. 



Finally, considering both the energy losses bi, 



bics we have 
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(17) 
(18) 



with the function yu(r) enclosing the whole spatial depen- 
dence. 

By substituting the above expressions into Ea. (|14[) we 
get the following equilibrium distribution for electrons 
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being /C = J U {-1}, bk = -afc/(fc + 1), if fc ^ -1, while 



C. Synchrotron spectrum 

The synchrotron spectrum of an electron gyrating in 
a magnetic field has its prominent peak at the resonance 
frequency 



v = 1^0 
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FIG. 4: Sky map at the frequency of 1 GHz for a realization 
of clumps distribution. For each clump, the circle radius is 
proportional to the logarithm of radio flux. 



with = 3.7 X 10^ Hz. This implies that, in practice, a 
i5-approximation around the peaks works extremely well. 
Using this frequency peak approximation, the synchrotron 
emissivity can be defined as 

M^,^ = ^(£;eH,r)^§^6,y„(£;,H,r). (21) 
d-bje dz^ 

This quantity is then integrated along the line of sight for 
the various cases to get the final synchrotron flux across 
the sky. 
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FIG. 5: DM synchrotron profile for the halo and unresolved 
substructures and their sum at 1 GHz for — 100 GeV 
and (aAv) = 3 x 10~^^cm^s~^. The astrophysical observed 
emission at the same frequency is also shown. The gray band 
indicates the angular region within which the DM signal from 
the host halo dominates over the signal from substructures 
modeled as in section III Al 



with del the distance between the observer and the clump. 
This can be rewritten as 
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(24) 



1. Single clump signal 



According to the description of previous section III Al 
let us consider a clump of mass md, whose center of 
mass is placed at i?ci and with a sufficiently small size. 
In this case it is possible to neglect the spatial variation 
of ^{f) inside the clump itself, and thus the flux can 
be calculated as: 



Iy{v,Rc\) = 



1 



drj,{v,R,i + r) , (22) 



clump 


del 


rc\ 


Pel 


B 


flux density at 1 GHz 


# 


kpc 


kpc 


GeVc-^cm-^ 


AiG 


GeVcm-^s-'Hz-i 


1 


14.2 


0.180 


6.51 


0.0962 


1.70 X 10"^^ 


2 


4.71 


0.181 


6.50 


0.320 


4.55 X 10"^^ 


3 


5.50 


0.188 


6.404 


3.08 


2.66 X 10"^^ 



TABLE H: Parameters of the example clumps chosen in Fig. 



El 



with 



GeV cm-2s-i Hz^i 



= 2.57 X 10 



kpc 



-12 



100 GeV 
kpc 



10-26, 
Pel 



(25) 



Fig. [3] shows some examples of signal, produced by 
three clumps of our simulation. An important feature 
to notice is that the synchrotron signal sensibly depends 
on the magnetic field both in the normalization and in 
the covered frequency range. In particular, the signal 
frequency cutoff, remnant of the energy spectrum cutoff 
near m^, depends on B following Ea l20l 

Fig. |4]shows instead the positions and radio intensities 
for a realization of the clumps distribution with masses 
TTici > IO^Mq. It can be seen that all the clumps with 
a non negligible signal lie near the galactic plane where 
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most of the galactic magnetic field is concentrated. Few 
clumps are visible at high latitude just because of pro- 
jection effects, being located very near and slightly up 
or below the solar position with respect to the galactic 
plane. 



2. Diffuse signals 

The diffuse halo signal is similarly given by the integral 
along the line of sight of Eq. ([2T|) 



dldb 



cos b 
in 



(26) 



where (/, b) are coordinates on the sphere and s the line 
of sight coordinate. To calculate the total contribution 



from the substructures, instead, we have to sum over all 
haloes 



d^/.'J 



dl db 



cos b I dnici I ds s 



dUc 
dnir 



■(mei,r)/r(^^,r), 



with given by Eq.dlSl) and r = f{s, I, b). 



(27) 



Interestingly, the sum of the two diffuse contributions 
can be rewritten as 



2 7DM 



dl db 



cos b 
An 



(28) 



where 



•DM 



1 / m^c 
4 



3 -1 

cm-'s 



GcV 

y. ^J.{r)^Ak{m^) 



Ph/GeVc ^cm 



(r/rh)(l + r/rh)2 



Pcl/GcVc ^cm ^ 
(r/rh)(l + r/r^Y 



l-fe/2 



(29) 



Thus, from the point of view of DM annihilation the 
unresolved clumps signal behaves like a further smooth 
NFW component with the same scale radius of the halo 
profile, but with a different effective density pch = 



-3 



0.604 GeVc cm , and with an emissivity simply pro- 
portional to the density profile instead of its square. 

We see that the halo component dominates in the cen- 
tral region of the galaxy, where 



r 



< 



PCh 



r < 4.39 kpc 



(30) 



which corresponds to a disk of radius 27.3 degrees (see 
fig.©. 



IV. DM ANNIHILATION CONSTRAINTS 

The pattern and intensity of the DM radio map result- 
ing from the sum of the contributions from the smooth 
halo and unresolved clumps is shown in Fig. [6] for 
= 100 GeV and (ctaw) = 3 x lO-^^cm^s^^ Similar 
maps are obtained at different frequencies and different 
and {gav) to obtain DM exclusion plots. For our 
analysis we use a small mask covering a 15° x 15° region 
around the galactic center where energy loss processes 
other than synchrotron and ICS start possibly to be rele- 
vant. We include the galactic plane although this region 



has basically no influence for the constraints on the DM 
signal. 

In FiglT] we show the radio constraints on the DM an- 
nihilation signal in the m^-{aAv) plane for various fre- 
quencies and various choices of the foreground. Several 
comments are in order. First, we can see that, as ex- 
pected, the use of the haze at 23 GHz gives about one 
order of magnitude better constraints with respect to the 
synchrotron foregrounds at the same frequency. How- 
ever, using also the information at other frequencies al- 
most the same constraints can be achieved. This infor- 
mation in particular is complementary giving better con- 
straints at lower DM masses. This is easily understood 
since a smaller DM mass increases the annihilation sig- 
nal (oc w~^) at smaller energies, and thus smaller syn- 
chrotron frequencies. In particular, the constraints im- 
prove of about one order of magnitude at ~ 100 GeV 
from 23 GHz to 1 GHz while only a modest improvement 
is achieved considering further lower frequencies as 0.1 
GHz. This saturation of the constraints is due to the fre- 
quency dependence of the DM signal, that below 1 GHz 
becomes flatter than the astrophysical backgrounds so 
that the fraction of contribution from DM is maximal at 
about 1 GHz. Further, the constraints show a threshold 
behavior given basically by Eql^D] which settles a maxi- 
mum emitted radio frequency for a given DM mass m^. 
This threshold behavior is, for example, clearly seen at 
23 GHz in the right panel of FiglT] where only for masses 
above 40 GeV the cross section is constrained. 
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FIG. 6: Sky map of the galactic radio signal generated by 
the DM smooth halo and unresolved clumps at the frequency 
of IGHz for = 100 GeV and {aAv) = 3 x 10"2®cm^s-\ 
The peculiar shape of the signal, pinched approximately at 
±30° and ±60°, reflects basically the structure of the mag- 
netic field as seen in projection from the Solar System, where 
the observer is located (compare also with Fig[2} . The galac- 
tic center region and the first few magnetic arms are visible as 
regions of high magnetic field intensity and hence prominent 
synchrotron emission. 



Notice that although the astrophysical background 
which we compare with at 1 GHz is an interpolation, 
the derived constraints are still valid given the smooth 
behavior and the broad frequency extent of the DM sig- 
nal, which does not exhibit narrow peaks at particular 
frequencies. However, effective measurements have been 
performed for example at 408 MHz and 1.4 GHz (see 
[40j). Quoting our constraints at these exact frequencies 
would change the results only slightly. 

The DM signal has thus a broad frequency extent and 
also below 1 GHz is still relevant. This is a potential 
problem for the DM interpretation of the WMAP Haze 
given that, in the Haze extraction procedure, the ob- 
served radio emission at 408 MHz is used as template 
of the synchrotron background. In fact, naively, a DM 
signal at 23 GHz should be relevant at 408 MHz as well, 
unless either the DM mass or the magnetic field is so 
high to shift the DM contribution to higher frequencies 
and making it negligible at 408 MHz. 

The second relevant point to notice is that the con- 
straint depends quite sensibly on the magnetic field as- 
sumptions. The constraints we obtain with the TT model 
are generally almost two orders of mag nitude weaker with 
respect to the results reported in [2l|. They are instead 
more in agreement with |22| where Galprop has been em- 
ployed to calculate the DM synchrotron signal. For a 
closer comparison with [2l[ we choose, as they do, a con- 
stant magnetic field of 10 ^G although still keeping the 
Galprop ISRF model. Even in this case our derived con- 
straints are a factor of 5 weaker (despite the inclusion 
of the contribution from substructures). The remaining 
factor of 5 can be finally recovered usinga constant ISRF 
with C/rad = 5eV/cm^ as assumed in [2l|. In this case. 



in fact, the smaller values of /7rad reduces the ICS losses 
enhancing in turn the synchrotron signal. It should be 
said however that, while the magnetic field normaliza- 
tion is still quite uncertain, the ISRF is instead more 
constrained and a large variation with respect to the Gal- 
prop model seems unlikely. 

The constraints shown in Fig|7] extend down to 
~ 10 GeV, which is somewhat the mass limit for a con- 
servative analysis. It is clear that for low masses the con- 
straints come more and more from lower frequencies: For 
example for a WIMP of 30 GeV the data at 100 MHz are 
2 orders of magnitude more constraining than the data 
at 10 GHz. However, extremely low frequencies are not 
experimentally accessible. For a WIMP of 1 GeV, from 
Eql^D] with a magnetic field of 0{fj.G) only frequencies 
< 10 MHz would be useful to place constraints on the 
DM signal. Although observations at this frequency exist 
[40| . in general the survey sky coverage is quite incom- 
plete and the data quality is non-optimal. Observations 
in this very low frequency range should substantially im- 
prove with the next generation radio arrays LOFAR and 
SKA. WIMP masses below 1 GeV still would produce 
observable synchrotron radiation at the galactic center 
where the magnetic field is likely much higher than ^G 
scale (possibly O(mG) ). This kind of analysis would be 
however quite model dependent and would face further 
background uncertainties. 



V. SINGLE CLUMPS DETECTABILITY 

To have a reliable estimate of the sensitivity to a single 
clump detection diffusion effects cannot be neglected. Al- 
though the integrated synchrotron clump signal is given 
by Ea. (P5)) . the clumps appear extended rather than 
pointlike with a dimension typically of several degrees. 
As a reasonable approximation we can assume that the 
signal is spread over an area of radius equal to the diffu- 
sion length of the electrons / £> = y^K{Ee)Tioss{Ee), where 
K is the diffusion coefficient and rioss is the energy loss 
time given by Eq. fTT)) . We use for K the Galprop model 
M 



K = Kn 



EeO 



(31) 



with a reference energy Eeo = 3 GeV, a Kolmogorov 
spectrum 5 = 1/3 and Kq = lO^^cm^/s. 

Taking as reference the parameters of a very bright 
clump like the ^^3 in table HIl we get (Ipc = 3 x lO^^cm) 



ID 



Kx) t" w 1 kpc , 



(32) 



for i^e ~ 10 GeV and for a radiation density Uiad ~ 
5eV/cm'^. The energy losses are basically dominated by 
ICS thus the result is almost independent of the magnetic 
field value. Moreover the dependence on the electron en- 
ergy and the radiation density itself is very weak. Of 
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FIG. 7: (Left) Constraints in the m^-{(JAv) plane for various frequencies without assuming synchrotron foreground removal. 
(Right) Constraints from the WMAP 23 GHz foreground map and 23 GHz foreground-cleaned residual map (the WMAP Haze) 
for the TT model of magnetic field (filled regions) and for a uniform 10 ^G field (dashed lines). 



course the clumps will have a certain profile peaked in 
the center and will not be perfectly smoothed all over 
Id . However the dilution of the signal in the much larger 
volume with respect to the region of emission makes it 
quite hard to detect the clump. We can consider for 
example the signal from a very bright clump at a dis- 
tance of 5 kpc with a flux of lO^^^GeVcm^^s^^Hz^^ 
at 20 GHz, corresponding approximately to the char- 
acteristics of clump #3 in FiglS] With a dilution over 
1 kpc, the clump emission is seen under a steradian 
A = -nc? ~ 0.1 sr with a « 10°, giving a diffuse clump 
flux of lO-^iQeVcm-^s-iHz-isr-i. The WMAP sensi- 
tivity of about 10 /zK translates into a flux sensitivity^ 
of ~ 10~^^GeVcm~^s~^Hz~^sr~^, meaning that the ex- 
pected, optimistic signal is about 3 order of magnitude 
below the reach of the current sensitivity. The situa- 
tion is only slightly better at 150 MHz where the ex- 
pected LOFAR sensitivity is 50 mK Q i.e. ^ 2 x 
lO-i^GeVcm-^s-iHz-isr-i. 



The chance of clump radio detection seems thus quite 
poor even with the next generation experiments. On the 
other side, the fact that the signal is anyway extended 
and not pointlike makes the clump signal not really com- 
plementary to the diffuse component sharing the same 
systematics with a much fainter signal. It is likely thus 
that the a role for DM investigations in the radio will be 
played basically by the diffuse signal. 



^ At radio frequencies the Rayleigh- Jeans law Fv = j <?k^T is 
employed to translate fluxes into brightness temperatures 



VI. SUMMARY AND CONCLUSIONS 

Using conservative assumptions for the DM distribu- 
tion in our galaxy we derive the expected secondary radi- 
ation due to synchrotron emission from high energy elec- 
trons produced in DM annihilation. The signal from sin- 
gle bright clumps offers only poor sensitivities because of 
diffusion effects which spread the electrons over large ar- 
eas diluting the radio signal. The diffuse signal from the 
halo and the unresolved clumps is instead relevant and 
can be compared to the radio astrophysical background 
to derive constraints on the DM mass and annihilation 
cross section. 

Constraints in the radio band, in particular, are com- 
plementary to similar (less stringent but less model de- 
pendent) constraints in the X-ray/gamma band [ssl . [56j 
and from neutrinos (57| . Radio data, in particular, are 
more sensitive in the GeV-TeV region while neutrinos 
provide more stringent bounds for very high DM masses 
( > 10 TeV). Gammas, instead, are more constraining 
for < 1 GcV. The combination of the various obser- 
vations provides thus interesting constraints over a wide 
range of masses pushing the allowed window significantly 
near the thermal relic possibility. 

More into details, we obtain conservative constraints at 
the level of {(Jav) ^ 10"^'^ cm'^s"^ for a DM mass = 
100 GcV from the WMAP Haze at 23 GHz. However, de- 
pending on the astrophysical uncertainties, in particular 
on the assumption on the galactic magnetic field model, 
constraints as strong as {<jav) ~ 10^^^ cm'^s^^ can be 
achieved. Complementary to other works which employ 
the WMAP Haze at 23 GHz, we also use the informa- 
tion in a wide frequency band in the range 100 MHz- 100 
GHz. Adding this information the constraints become 
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of the order of {cr^v) ~ 10~^^cni^s~^ for a DM mass 
rriy. = 100 GeV. The multi-frequency approach thus gives 
comparable constraints with respect to the WMAP Haze 
only, or generally better for ^100 GeV where the 
best sensitivity is achieved at ~ GHz frequencies. 

The derived constraints are quite conservative because 
no attempt to model the astrophysical background is 
made differently from the case of the WMAP Haze. In- 
deed, the Haze residual map itself should be interpreted 
with some caution, given that the significance of the fea- 
ture is at the moment still debated and complementary 
analyses from different groups (as the WMAP one) miss 
in finding a clear evidence of the feature. In this respect, 
the multifrequency approach will be definitely necessary 
to clarify the nature of controversial DM signals as in the 
case of the WMAP Haze. Progresses arc expected with 
the forthcoming data at high frequencies from Planck and 
at low frequencies from LOFAR and, in a more distant 
future, from SKA. These surveys will help in disentan- 
gling the various astrophysical contributions thus assess- 



ing the real significance of the Haze feature. Further, the 
low frequency data in particular, will help to improve 
our knowledge of the galactic magnetic field. Progresses 
in these fields will provide a major improvement for the 
interpretation of the DM-radio connection. 
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